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Abstract Enormous progress has been made on observing stellar magnetism 
in stars from the main sequence (particularly thanks to the MiMeS, MAGORI 
and BOB surveys) through to compact objects. Recent data have thrown into 
sharper relief the vexed question of the origin of stellar magnetic fields, which 
remains one of the main unanswered questions in astrophysics. In this chapter 
we review recent work in this area of research. In particular, we look at the 
fossil field hypothesis which links magnetism in compact stars to magnetism in 
main sequence and pre-main sequence stars and we consider why its feasibility 
has now been questioned particularly in the context of highly magnetic white 
dwarfs. We also review the fossil versus dynamo debate in the context of 
neutron stars and the roles played by key physical processes such as buoyancy, 
helicity, and superfluid turbulence, in the generation and stability of neutron 
star fields. 

Independent information on the internal magnetic field of neutron stars 
will come from future gravitational wave detections. Coherent searches for the 
Crab pulsar with the Laser Interferometer Cravitational Wave Observatory 
(LICO) have already constrained its gravitational wave luminosity to be < 2 % 
of the observed spin-down luminosity, thus placing a limit of < 10^® C on 
the internal held. Indirect spin-down limits inferred from recycled pulsars also 
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yield interesting gravitational-wave-related constraints. Thus we may be at the 
dawn of a new era of exciting discoveries in compact star magnetism driven 
by the opening of a new, non-electromagnetic observational window. 

We also review recent advances in the theory and computation of mag¬ 
netohydrodynamic turbulence as it applies to stellar magnetism and dynamo 
theory. These advances offer insight into the action of stellar dynamos as well 
as processes which control the diffusive magnetic flux transport in stars. 

Keywords Magnetic helds • Main sequence stars • white dwarfs • neutron 
stars 


1 Introduction 


It was Hale 


one of the greatest astrophysicist of the twentieth century, 
who built the first spectroheliograph and used it to establish that sunspots are 
magnetic and grouped in pairs of opposite polarities. Furthe r pion eering work 
on solar magnetism conducted by Hale and collaborators in 19191 revealed an 
East-West direction of polarity in the sunspots’ magnetic fields exhibiting a 
mirror symmetry with respect to the solar equator. Such polarity was observed 
to undergo inversion according to the 11 years solar cycle. This is commonly 
known as the “Hale-Nicholson law”. 

The first detectio n of a magnetic field in a sta r other than our own Sun, 78 
Vir, was obtained in 19471 by Babcock. In 1958 , Babcock published the first 
catalogue of magnetic stars and in 1960 he discovered what is still today the 
most magnetic main sequence star known, HD 215441 (~ 3.4 x 10^ G). This 
became known as “Babcock’s star”. Surface fields discovered in more recent 
years that rival i n strength that of Babcock’s star are those of H D 154708 


2 . 45 X lO^G: iHuhris- et 3.1.11200511 . HD 137509 2.9 x lO^G 


19951: lKochukhoyll2006ll and HD 75049 (- 
Elkin et al. 12010 1. 


3 X 10^ G; Frevhammer et al. 


Mathvj 

" 20081 : 


Following the detections of strong fields in main sequence stars, BlackettI 


(1943) suggested that if magnetic flux is conserved, some white dwarfs should 
exhibit magnetic fields of up to 10^ — 10® G. However, spectroscopic surveys 
aimed at dete ct ing magnetic fields in white dwarfs yielded negative results 
(|Prestonlll970ll . |k emd (jl970ll argued that electrons spiralling in a magnetic 
field would emit linearly and circularly polarised radiation that should be 
detectable in the continuum of strongly magnetic white dwarfs. A spectropo- 
larimetric survey of white dwarfs led to the discovery of strong circular po¬ 


larisa tion in the continuum of the white dwarf Grw-|-70°8247 (|Kemp et al 
1970lb 


Baade and Zwickv ( 1934h first suggested that some stars could be made up 


of neutrons and that a supernova could be the resul t of a rapid trans ition of 
a normal star into a neutron star. iBaadel ( 1942h and Minkowski ( 1942 1 found 
unusual emission arising from the central parts of t he Grab Nebula. Later 
on radio pulsations in this nebula were discovered by IStaelin and Reifenstein 
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(1968). Such radio emission had been predicted bv IShklovskv ( 1958 i as caused 
by relativis t ic elec trons spiralling around magnetic field lines. 

IWoltier ( 1964h first proposed that under magnetic flux conservation, if a 
star contracts to the density expected in a neutron star, then the surface field 
strength could be amplihed to values of up to 10^ "*^ — 10^^ G. T he first pulsar, a 


highly magnetised rapidl y spinning neutron star (|Pacinilll968f) . was discovered 


in 1967 by Jocelyn Bell (jHewish et al.lll968f) . 


In this chapter we review progress made on the origin of magnetic fields 
in stars. The origin of magnetic fields is still a major unresolved problem in 
astrophysics. The “fossil field” hypothesis is often invoked to link magnetism in 
compact stars to magnetism on the main sequence (see Sect. 13. H and Sect. 14.1|) . 
In the fossil scenario, some fraction of the magnetic flux of the progenitor 
star is conserved during the collapse process, because the stellar plasma is 
highly conducti ng and hence , by Faraday’s law, the magnetic held is ‘frozen 
in’ to the fluid (IW oltieil r964 iRnderman 1972 : Braithwaite and Soruitl 2nn4at 

' rhel2 


iFerrario and Wickramasingha200~ . 20071 2008a b Under these circumstances, 

the field strength B scales with radius R of the star as B (x R~^. A main 
sequence B star with radius 5Rq and field ~ 1,000G compresses to give a 
neutron star with a field of ~ 10^^ G; a main sequence A star with radius 
1 . 6 i ?0 and field ~ 1 , 000 G compresses to give a magnetic white dwarf with 
a field of ^ 10^ G. The fossil field scenario is indeed quite attractive and can 
explain the existence of magnetic fields even in the strongest magnets in the 
universe, the so-called “magnetars”. However, the feasibility of the fossil field 
hypothesis does have its problems and has been recently questioned, in the 
context of highly magnetic white dwarfs, on the basis of recent observational 
results related to their binarity (see Sect. 13.21) . The fossil versus dynamo¬ 
generated fields debate for neutron stars is analysed in Sect. 14.1.11 Alternative 
explanations for the origin of fields in stars are presented in Sect. 13.21 and 


In this chapter we also review recent progress in magnetohydrodynamics 
(MHD) turbulence theory and computation as it applies to stellar magnetism 
(see Sect. 0 , which underpins the operation of stellar dynamos, and controls 
diffusive magnetic flux transport in stars. 


2 Magnetism in non-degenerate stars 

Magnetic fields in main sequence stars have been measured mainly via spec- 
tropolarimetry. Different techniques measure different held properties, e.g., 
linc-of-sight component, volume-averaged held strength, or dipole moment. 
A thorough review on magnetic held measurements of non-degenerate stars 
across the Hertzprung-Russell diagram can be found in this book in the chap¬ 
ter by Linsky & Scholler. 

Direct measurements of magnetic helds in the chemically peculiar main 
sequence Ap and Bp stars, which form about 10% of stars in the 1.6 — SMq 
range, have revealed the existence of strong (^ 3 x 10^ — 3 x 10^ G) large scale 








































4 


Ferrario, Melatos, Zrake 


fields (e .g. Auriere et, al. 2007 : Donati and Landstreetll200^ . Landstreet et al. 
(|2007all found that the majority of magnetic objects have an average longitu¬ 
dinal field (line of sight component of the surface field) B > 0-25 kG in every 
1 Mpt m ass bin between 2Mq and 9Mq ( Baenulo et alJlioOfia : Landstreet et alJ 


2007all . The lack of fields below 300 G, value that corresponds to the strength 


at which the magnetic field is in equipartition with the ga s pressure in the stel- 
lar photosphere, is not a detection threshold effect since Auriere et al.l ( 2010ll 
have set a 3 ct upper limit of longitudinal fields dowrr to 1 — lOG depend¬ 
ing on stellar brightness. The absence of magnetic stars below this cutoff has 
been referred to as “the Ap/Bp magnetic desert”. This indicates that Ap and 
Bp stars are not simply the high field tail of a continuous field distribution. 

also find that the field strengths seem to show some 


Landstreet et al. 


decline with stellar age but the field iircidence does not. Such large scale fields 
are observed throughout the main sequence phase and, in more recent years, 
have also been obse rved in a small number of stars on the red giant branch 
( Auriere et ahllioO^ . 


The recent detections of sub-gauss fie lds in Vega ( Lignieres et al. 2009ll . 
and p ossibly in the Am stars Sirius A l)Petit et al.l 1201111 and /3UMa and 
9 Leo ( Blazere et al. 2014h have unveiled a new class of magnetic stars which 
is at the 1 gauss end of the magnetic desert, thus potentially suggesting a di¬ 
chotomy bet ween strong and ultra -weak magnetic fields among intermediate- 
mass stars ( Li gnieres et al.l 20 14ll. This in triguing observational result has 
promoted iBraithwaite and Cantiell^ ( 2013 1 to investigate the origin of these 
low fields. They argued that such fields could be the remnants of fields already 
present or formed during or immediately after the star formation stage. Hence, 
these fields would still be evolving on a timescale that is comparable to the 
age of the star. According to these studies, all intermediate and high mass 
stars lacking strong fields should display sub-gauss field strengths that would 
slowly decline over their main sequence lifetime. There has been some r ecent 
effort to test this pre diction bv lWade et al.l (2014a), Neiner et al.l ( 2014bll and 


Neiner et al.l (|2014afl but without any clear conclusion. 


lAlecian et al.l ( 2008ll conduct ed a spectropolarimetric study of Herbig Ae/Be 
(HAeBe) objects (Herbig 1960l) . HAeBe stars are pre-main sequence stars 
of 2 — 15 Mq which are still embedded in their protostellar gas-dust enve- 
lope and exhibit em i ssion lines in spectra of type A/B. The o bservations of 


lope ana exnibit ern i ssion lines m spectra ol type A/lj. ine o bservations ol 
lAlecian et al. (l2008l) . Alecian et al.l (12009 ). Hubrig et al. ( 20091) and more re- 


cently Alecian et al.l (2013a) and Hubrig et al.l ( 2013 ) have revealed that about 


7% of HAeBe objects are magnetic. These studies have also found that they 
display large scale dipolar fields of strength comparable to that of Ap and Bp 
stars, under the assumption of conservation of magnetic flux. Interestingly, 
although some magnetic HAeBe stars are in the totally radiative evolution¬ 
ary phase while others have already developed convective cores, they seem 
to share the same magnetic field strength and structure thus indicating that 
the convective core i s not responsible for generation or destruction of fields 
( Alecian et al. 2013b ). 
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The M agnetism in Massive Stars (MiMeS.IWade et alJl201lli project, the 


ag 

MAGORI (jHubrig et alJ 120111 : ISchoUer_et_ajjj20 


in OB star (BOB) Collaboration (|Hnbrig et al 


ID p roject and the B fields 


2ni4D have conducted large 


surveys of bright and massive Galactic stars of spectral types B and O. It is 
now clear that about 10% of all stars with radiative envelopes in the mass range 
1.5— 50 M© possess large scale mostly dipolar magnetic fields (IGrunhut and Wade 
201, iD . Even more interestingly, the MiMes project has revealed that all of the 
basic field characteristics do not vary significantly from the coolest spectral 
types FO 1.5 M©) to the hottest spectral types 04 50 M©) stars thus 

support ing a common form ative scenario over a very large range of stellar 
masses ( Wade et al.ll2014b l. Furthermore, the upper dipolar field limits placed 
on the undetected magnetic O stars sample studied by the MiMeS collabora- 


tion are 40 G at 50% confidence, and 130 G at 80% confidence (Wade and the MiMeS Gollaboration 


2ni4D . This result seems to indicate that the field distribution of massive stars 
may also have a magnetic desert, simila r to that observed i n intermediate- 
mass stars. However, the recent work of iFossati et al. ( 2014 1. who detected 
weak longitudinal fields of < 30 G in the early B-type main sequence stars 
/? CMa and e GMa, appears to support a more continuous distribution of fields 
in massive stars. They also claim that weak fields in massive stars could be 
more widespread than currently observed because of the numerous observa¬ 
tional biases associated with the detection of weak fields in massive stars using 
current instrumentation and techniques. 


The vexed question on the origin of fields in main sequence stars is still 
unanswered. The two main hypothesis are the fossil field, according to which 
magnetis m in stars is a relic of the inte rstellar field from which the star was 
born (e.g. Borra et al. 19821 IMossI[ 200lh . and dynamo action taking place in 
the rotating cores of main sequence stars. Both theories are unable to ex¬ 
plain why only a small subset of main sequence stars is magnetic, although 
the fossil field theory does assume that differences such as initial density and 
magnetic field strength in interstellar clouds could be responsible for this. How¬ 
ever, it is peculiar that observations have never revealed the presence of mag¬ 
netism in both components of a binary system. The Bina MIcS (Binarity and 
Magnetic Interactions in various classes of Stars) project ( Alecian et al.l 2014i 
Neiner and Aleci^ 201,3D is a collaborative endeavour that has been recently 
setup to investigate the phenomenon of magnetism in close (orbital periods 
shorter than 20 d) binary systems. The study of the magnetic properties of 
these systems, at a detection limit of 150 G, will give us crucial information 
on the origin of fields i n magnetic main se quence stars which are well known 
to be rare in binaries ( Carrier et al.l [2002 1. Since the components of binaries 
share the same history, the study of the two stars will help us unravel the 
importance of their physical parameters and their studies should enable us to 
distinguish effects caused by initial conditions at formation from those caused 
later on by evolutionary processes. Preliminary studies of 314 massive stars 
in close binaries have revealed the presence of fields in 6 systems yielding an 
incidence of magnetism of < 2% which is up to 5 times lower than what is ob- 
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( 2002 ^ which 


served in single stars. This conhrms the results of iCarrier et al 
were based on a smaller sample of A-type stars. 

In order to explain the rarity of short period binaries containing magnetic 
main sequence stars and also the fact that magnetism is only pr esent in a 
small percentage of single main sequence stars, Ferrario et al. ( 2009ll proposed 
that magnetic fields in main sequence stars could form when two proto-stellar 
objects merge late on their approach to the main sequence and when at least 
one of them has already acquired a radiative envelope. T he N-body simulations 
of pre-main sequence evolution of iBailton et al. ( 20l4l have indicated that in 
high initial cluster densities the number of collisions between stars is twice as 
high as that of stars already on the main sequence. Furthermore, they find that 
massive stars generally form through the merging of lower mass stars. Thus, 
the expectation is that the incidence of magnetism sho uld increase with m ass, 
which seems to be demonstrated by the observations of Power et al. for 

a volume limited sample of A and B type stars with M < 4Mfri . Observations 


also s eem to be in general agreement with the predictions of iFerrario et al 


( 2 OO 9 II that the Ap and Bp stars should have pre-main sequence progenitors 
with similar magnetic field flux and structure. _ 

Another merger scenario has been proposed bv lTutukov and Fedoroval ( 2010ll . 
They propose that the coalescence of the two components of a close binary 
system with masses in the range 0.7 — 1.5 Mq, which are expected to have 
convective envelopes and strong magnetic fields, would be responsible for the 
formation of Ap and Bp stars with M < 3Mq. In this picture, the magnetic 
Ap and Bp stars with radiative envelop es owe their strong magnetic fields to 
progenitors with convective envelopes. Tutukov and Fedorova ( 2010ll specu¬ 
late that systems such as the W UMa-type contact binaries could represent 
the precursory phase of these merging events. 

However all of these predictions are so far only of a very qualitative nature 
and need to be supported by detailed quantitative calculations in the future. 


3 Magnetic white dwarfs 


Magnetism in white dwarfs is revealed either through Zeeman splitting and cir¬ 
cular polarisation of spectral lines or, at very high field strengths, as continuum 
polarisation. The magnetic white dwarfs, which represent about 8 — 13 % of 
the total white dwarf population ( Liebert et al.ll2003l: iKawka et al. 2007^ and 
exhibit polar field strengths of ^ 10^ — 10® G, have been discovered mainly in 
optical sky surveys. At spectral resolutions of about 10 A, the Zeeman triplet 
in the lower members of the Balmer and Lyman series can be easily recog¬ 
nised at fields of ~ 10® — 10^ G when the splitting in intensity spectra is 
smaller than other broadening factors, such as pressure broadening. At higher 
spectral resolutions, it is possible to detect fields down to ^ 10® G. However, 
circular spectropolarimetry is a much more sensitive observational technique 
that is used to measure even lower (< 10® G) fields, since it can independently 
measure the and a~ oppositely polarised components of the Zeeman split 
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absorption features. On the other hand, due to the faintness of white dwarfs, 
fields below ^ 3 x 10"^ G have onl y been recently detected . Spect r opolarimetric 


observ ations on the ESO VLT bv lAznar Cnadrado et al. ( 20041 1 : Ijordan et ^ 
(1200711 revealed fields down to ~ 10^ G in about 10% of the surveyed white 


dwarfs iLandstreet et al.l (|2012ll have confirmed this finding by conducting a 
survey of objectsrandomly drawn from a list of nearby cool {Tes < 14, 000 K) 
WDs. 

Observations seem to indicate that there is a p aucity of white dwarfs wit h 
fields in the intermediate field range 10® — 10® G ( Kawka and Venneal2012|l . 
reminiscent of the magnetic desert of Ap and Bp stars. However, this finding 
has not been fully corroborated by observations and future surveys may find 
objects also in this field range. In this context we need to stress that while 
magnetic helds in bright main sequence and pre main sequence stars have been 
mainly found via circular polarisation observations, only a small percentage 
of magnetic white dwarfs has been discovered using this method. This non- 
systematic approach has created, over the past 30 years, a sample of stars 
whose properties are difficult to investigate since the observational biases are 
not fully understood and thus are difficult to remove. 

Another interesting characteristic of strongly magnetic white dwarfs is that 
their mean mass (O.85±O.O4M0) is higher than that of non-magnetic or weakly 
magnetic white dwarfs (0.593 ± 0.002 Mq, see the chapter on magnetic white 
dwarfs in this book). This indicate that their progenitors are more massive 
than t hose of ordinary white dwarfs. Interestingly, the observations of lVenn^ 


(Il999ll of the EUVE/Soft X- ray selected ultra-massive white dwarfs (M > 
1.2 Mq) have found that ~ 20 % of them are strongly magnetic. Thus it appears 
that the incidence of magnetism in the high field group increases with white 
dwarf mass and hence with progenitor mass, unless they result from a merger 
(e.g. EUVE ,10317-853 Barstow et al. fl995l: iFerrario et al. 1997 1. 


3.1 Origin of helds in highly magnetic white dwarfs: the fossil hypothesis 

According to the “strong” fossil held hypothesis, magnetic stars originate in 
gas clouds whose helds are at the upper end of the held d istribution observed in 
the interstellar medium (10“® and 10“^ G, HeilesI 1997h. One of th e very hrst 
articles on star formation and hux freezing is tha t of Mestell l|l96fill . Briehy, in 
the simple model proposed bv lTout et al. ( 2004 1. a star of mass M that col¬ 
lapses from the interstellar medium entraps a magnetic hux oc Mx 

where Hism is the interstellar (primordial) held. The variations expected in the 
interstellar magnetic held will determine the distribution of magnetic huxes 
in the protostars. The magnetic hux lines would freeze in the radiative upper 
layers of the emerging star that will then evolve towards the main sequence. 
Assuming perfect magnetic hux conservation, a contraction by a 10^ factor 
could in principle create helds of the order of 10® G in a main sequence star. 
However during the cloud’s collapse most of the initial magnetic hux would 
be lost via ohmic dissipation or ambipolar diffusion, so that much lower helds 
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than this upper limit can be realised. The survival of significant large-scale 
magnetic flu x through the pre-main sequence evolution has been addressed by 
Mos j (l2nn.ih . 

According to the “weak” fossil field hypothesis, stellar fields could just be 
the remnants of those generated by the dynamos of active pre-main sequence 
stars or by late mergers of proto stars with at least one of the two having a 
radiative envelope (see Sect.and iFerrario et ahllioOflll . The magnetic field flux 
would then be conserved during the subsequent evolution from main sequence 
to the white dwarf stage. 

Under the fossil field hypothesis, the magnetic Ap and Bp stars in the mass 
range 1.5 —8 Mq would be the progenitors of the highly magnetic white dwarfs 
{B > 10® G). The paucity of white dwarfs with fields in the range 10® — 10® G 
could be interpreted as being related to the lower field cut off of Ap and Bp 
stars. 

The main issue is whether magnetic fields can survive complex phases of 
evolution when a star develops convective and ra diative zones that contract 
and expand in size with time. iTout et al.l (1200411 point out that a dynamo¬ 


generated field in the convective regions of a star is transient and has no 
large-scale component. A star can conserve its primordial fossil field as long as 
this star possesses a radiative zone throughout its life, because the field would 
be pumped out of any dev eloping convective region into an adjacent radiative 
region. The conclusions of Tout et ^ ( 2004ll is that a poloidal fossil field that 
is present in a main sequence star can appear in the compact star phase with 
a similar field structure and with its magnetic flux nearly conserved. 

There have been many semi-analytical and numerical calculations that have 

addressed the issues concerning the evolution of a m ag netic field during the _ 

main se quence phase. Braithwaite and Spruit ( 2004ljl and iBraithwaite and Nordlund 
(l2006all have shown that stable magnetic fields with roughly equal poloidal 
and toroidal field strengths can exist in the radiative interior of a star and 
their exterior appearance would be that of a dipole with minor contributions 
from higher multipoles (see also iMestel and MossI l201flll . Interestingly, they 
find that Ohmic dissipation diffuses the field outward over time so that the 
field strength at the surface of the star increases while the field structure in 
the stellar interior would switch from being mainly toroidal to poloidal. Hence 
they predict an increase in the surface field strength with the age of the star. 

In order to test these theoretical findings, it is necessary to study a sample of 
Ap and Bp stars with known age and fraction of the main sequence evolution 
completed. However , it is quite difficult to estimate age-related quantities for 
magnetic stars (see iBagnulo et al.l l2006bll . The evolution of magnetic fields 
from the zero-age main sequence (ZAMS) to the terminal-age main sequ ence 
(TAMS) has been invest igated observationally bv IBagnulo et al. ( 2006ljl and 
Landstreet et al.l ( 2007lJl through the study of Ap stars in clusters, since the 


age of a cluster can generally be established to better than about ±0.2 dex. 
They find that magnetic fields exist at the ZAMS phase for stellar masses 
2 — 5M0 when stars have fractional ages (that is, fraction of the main se¬ 
quence evolution completed) below about 0.05 and for fractional ages of less 
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than about < 0.10 for masses up to 9Mq. This is consistent with the exis¬ 
tence of fields in the late pre-main sequen c e stars Herbig AeBe (see section 
El). However, the study of iLandstreet et ahl (|2007lJl also reveals that the evo¬ 
lutionary time of magnetic fields is less straightforward to interpret. For stars 
with M > 3 Mq they find that the field strength dec reases on a timescale of 
about 23 X 10^ yr, in agreement with iKochukhovI ( 2006h . However it is not clear 
whether the magnetic flux really decreases or whether it is the conservation 
of magnetic flux that causes the surface field to decrease as the star expands. 
Interestingly, for stars with M <,3M0 they find no conclusive evidence of 
field strength or magnetic flux reduction on a time scale of a few 10® yr. 


Support for the fossil field hypothesis also come s from spectropolarimetric 


studies of fifty red giants by lAuriere et al.l (120131) . This work has revealed 


the existence of four magnetic giants that have been identified as probable 
descendants of Ap/Bp stars. 


Po pulation synthesis calculations carried out bv lWickramasinahe and Ferrario 
( 2 OO 5 I) have shown that starting with a distribution of magnetic fields on the 
main sequence, as observed for the Ap and Bp stars and under the assumption 
of magnetic flux conservation, the predicted magnetic field and mass distribu¬ 
tion in white dwarfs are in general agreement with observations, provided that 
a modified initial to final mass relation is employed for magnetic white dwarfs 
(their “scenario A”). That is, in order to fit the mass distribution, one needs 
to assume that mass loss mechanisms are partly inhibited by the presence of 
strong fields. However, such models can at most yield an incid ence of mag¬ 
netis m of ^ 5% as compared to the o b serve d ^ 8 — 12% (see also lKawka et al. 


20071) . IWickramasinghe and Ferrariol (|20051) have shown that it is possible to 


achieve a better agreement with observations if one assumes that in addition 
to the actually observed fields in Ap and Bp stars, about 40% of stars more 
massive than 4.5 Mq (early B-type stars) have fields of 10 — 100 G and evolve 
into high fields magnetic white dwarfs (see Figured]). This is their “scenario B” 
which would only be viable if magnetic B-type stars do not have a “magnetic 
desert”, whic h seem s to be, albeit not conclusively, supported by the work of 


Fossati et al.l (20l3). 


Larger and more sensitive spectropolarimetric surveys of white dwarfs may 
confirm the classification of white dwarfs within either the very low (< 10^ G) 
or the very high field group (> 10® G), hence forming a bimodal field distri¬ 
bution. This would point to two different channels for field formation. It is 
interesting to note that th e weak magnetic fields which have rece ntly been de¬ 
tected in Vega and Sirius ( Lignieres et al. 20091 : IPetit et aklfeoill) . would scale 
up to a surface field of a few ^ 10^ G at the white dwarf stage under magnetic 
flux conservation. Such weakly magnetised white dwarfs would belong to the 
low field component of the bimodal white dwarf field distribution. 
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Fig. 1 Population synthesis calculations of lWickramasinghe and Ferrari^ (1200511 (their sce¬ 
nario B) for the field distributions of magnetic white dwarfs based on the fossil field hypoth¬ 
esis compared with observations 


3.2 Origin of fields in highly magnetic white dwarfs: the stellar merger 
hypothesis 


We have plotted in Fig. [5] the ratio of poloidal magnetic flux to mass of the 
highest field magnetic white dwarfs and of the most magnetic main sequence 
stars. This diagram shows that the two groups of stars share similar char¬ 
acteristics indicating that their magnetic fields may have a common origin. 
However, the possible evolutionary link between the two populatio ns, and thus 


the vi abilit y of the fossil field hyp o thesis , has been questioned by iTout et al 


(|2008f) and IWickramasin ghe et al.l 112014). Th e main argument hinges on the 
fact, first highlighted bv Ibiebert et al.l ( 200511 . that there is not a single ex¬ 
ample of a high held magnetic white dwarf {B > a few 1 0^ G) with a n on- 
degenerate companion star (generally an M dwarf, but see Ferrarioll2012li in 
a non-interact ing binary. Fu r ther s earches on a much larger sample of objects 
conducted by Liebert et al.l ( 20151) have conhrmed the hypothesis that mag¬ 
netic held in white dwarfs and binarity (with M or K dwarfs) are independent 
at a 9 CT level. However, we cannot invoke some as yet unknown physical mech¬ 
anism that could inhibit the formation of a strong magnetic held in a white 
dwarf when a companion star is present, because such a mechanism would also 
prevent the formation of Magnetic Cataclysmic Variables (MCVs) which are 
interacting binaries consisting of a magnetic white dwarfs with a low-mass red 
dwarf companion. This curious lack of detached white dwarf - non-degenerate 
star systems indicates that there are no known progenitors of MCVs (for more 
details see also the chapter on magnetic white dwarfs in this book). 

An merger s c enario for the generation of helds has also been proposed by 


Nordhaus et al.l l|201l[) . Here a low-mass star would be tidally disrupted by 


its proto-white dwarf companion during common envelope evolution to form 
an accretion disc. In this disc seed helds would be amplihed through tur¬ 
bulence and shear and then advected on to the object that will become an 
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Fig. 2 The ratio of magnetic flux to mass (f)-o/M for the most mag netic main sequence stars 
(squares) and white dwarfs (circles) (IWickramasinghe et al.ll2014fl 


isolated highly magnetic white dwarf. Garcia-Berro et al. ( 2012ll conducted 


three-dimensional hydrodynamic simulations of merging double degenerates 
and argued that a hot and differentially rotating convective corona would 
form around the more massive star. Their population synthesis calculations of 
double degenerate mergers are in general agreement with observations. Sim¬ 
ilar population synthesis calculations conducted for a wide range of merging 
astrop hysical objects have also been carried out by iBogomazov and TutukovI 

(I2nn9ll . _ _ _ _ 

Following the work of iTout et al.l ( 2008 ), Wickramasinghe et al. ( 2014 ) pro¬ 
posed that the fields are generated by an 17 dynamo within the common en¬ 
velope of a binary system where a weak seed poloidal held would wind up by 
differential rotation. The closer the cores of the two stars are drawn before 
the envelope is ejected, the stronger the hnal held of the star emerging from 
common envelope will be. Therefore the strongest helds occur when the two 
stars merge, forming an isolated highly magnetic white dwarf. If the two stars 
do not coalesce but emerge from the com mon envelope wh en they are about 
to transfer mass, they become the MCVs ( Tout et ^ 2008ll . 

Observations indicate that highly magnetic stars are typihed by (see Fig. 

El) 


10 


- 8.5 


< 


^p/M 

Gcm^g- 


< 10 


- 6.5 


( 1 ) 


where <Pp = R^Bp is poloidal magnetic hux and M and R are the total mass 
and radius of the star respectively. 

If the dynamo mechanism generates a magnetic held from differential ro¬ 
tation Af], then we have 
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Fig. 3 The evolution of the magnetic field components Sp and B) following 

a stellar merger which produces a differential rotation equal to the break up spin of the 
merged object. The decay of differential rotation AH follows the right hand axis while the 
fields follow the left hand axis. Toroidal fi eld decays unless ? 7 p > arp' and this determines 
the final ratio of toroidal to poloidal field dWickramasinghe et al.ll2Q14|j 


Toroidal and poloidal fields are unstable on their own ( Braithwaite 200tt . 
The process that gives rise to the decay of toroidal field leads to the generation 
of poloidal field with one component stabilising the other and thus limiting 
field growth until a stable configuration is reached . The results of Fig [5] have 
been interpreted by Wickramasinghe et al.l ( 2014ll as follows. If ifp, and 


E = Ep+E^ are the poloidal, toroidal and total magnetic energies respectively 
and r] = E/U, rjp = Ep/U and 77 ^ = E^/U are, respectively, the ratios of 
poloidal and toroidal magnetic t o thermal energy U, then by scaling to the 
observed maximum poloidal flux, Wickramasinghe et al.l ( 2ni4h find that 


Vp = 


10 


-8 




A y cm^ g“ 


( 3 ) 


where <Pp is the ratio of magnetic flux to stellar mass. The observational data 
in Fig. [2] show that the maximum r]p is independent of the mass and type of 
star. In non-rotatin g stars, a stable poloidal-toroidal configuration must satisfy 
( Braithwaitel 200*^ 

a?7^ < ?7p < 0.8?7, (4) 

where a ~ 1 is a buoyancy factor for main sequence stars. The left hand 
side inequality is caused by the stabilising effect of a poloidal field on the 
Taylor instability in purely toroidal fields. A lower limit to the poloidal field 
is determined by the relative importance of magnetic to gravitational-thermal 
energy through buoyancy effects. The upper limit is due to the stability of 
poloidal fields which requires they be not significantly larger than the toroidal 
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field. Braithwaitel ( 2009ll argued that the same inequalities are also likely to 
hold for stable fields in rotating stars. 

Fig.O where ( = 1/77 oc B, shows how Bp and B^ evolve as AH decreases 
from its maximum at Af2 = f7crit- Their single model parameter is chosen 
to give the observed maximum r]p Ri 10“®, so that Cp ~ 10“"^. Toroidal field 
is initially generated by the winding up of the seed poloidal field through 
differential rotation. As soon as the toroidal field is large enough, the poloidal 
field starts to grow. As long as the differential rotation is sufficiently large, 
the toroidal field continues to grow up to a maximum value. Then it decreases 
until it reaches equilibrium with the poloidal field and aij'^ = rip. The magnetic 
torque extinguishes Z\ f 7 to y ield a final object that rotates as a solid body. 
Wickramasinghe et al.l ( 2014ll find that the dynamical timescale for a 2 Mq 


star is of the order of 40 min so this evolution is completed in about 3.7 yr which 
is a lot smaller than the corresponding K elvin-Helmholtz timesca l e (2.3 x 


10® yr for a 2Mq star). The conclusion of Wickramasinghe et al. ( 2014ll 


IS 


that the final poloidal field is stable and proportional to the initial quantity 
of differential rotation, but is independent of the size of small initial seed field 

and. _ 

More sophisticated a — dynamo models are described in iNordhaus et al 


( 2007ll . although these calculations are aimed at understanding the different 
problem related to the sha ping mechanisms in pla netary nebulae. 

In any case, the work of Nordhaus_et.j|l. ( 201ll) highlights that for the dy¬ 
namo mechanism of IWickramasinghe et al. ( 2014 1 to be viable, the transport 
of a strong field generated in the envelope to the proto-white dwarf’s surface 
would req uire a diffusion coefficie nt > 10^^ — 10^^gcm“® which is unphysi¬ 
cal. Hence Nordhaus et al.l ( 201lll propose an alternative scenario whereby a 
companion star, of sufficiently low mass to avoid a premature ejection of the 
envelope, would be tidally destroyed by the gravitational field of the primary 
star and would form a disc extending all the way to the surface of the proto¬ 
white dwarf. Accretion would then transport the strong fields that are formed 
in the disc towards the proto-white dwarf where they would become anchored 
to its surface. 


The population synthesis calculations carried out bv iBriggs et al.l (|2015ll 


have given further support to the stellar m erging hypothesis fo r the origin of 
fields in the highly magnetic white dwarfs. Briggs et al.l ( 2015ll find that the 
mass distribution and the fraction of stars that merge during a common en¬ 
velope phase are in good agreement with the observations of magnetic white 
dwarfs for a wide range of common envelope efficiency parameter a. Fig.|4] 
shows the contributions from the various merging pre-common envelope pro¬ 
genitor pairs for a = 0.1. They find a theoretical incidence of magnetism 
among white dwarfs, for a = 0 .1 — 0.3, of about 13-19%, which is consistent 
with observational results ('e.g. lKawka et al.ll2007ri . 

For the sake of completeness, we remark that similar population synthe¬ 
sis calculations, but related to the viability of the core-degenera.te sc enario 
for Type la supernovae, have been conducted bv lllkov and Soke^ ( 2013ll (and 
references therein). 
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Fig. 4 Theoretical mass distribution of magnetic white dwarfs for a = 0.1. The contribu¬ 
tions are separated according to their pre-CE progenitors. AGB= Asymptotic Giant Branch; 
MS= Main Sequence; RGB= Red Giant Branch; CS= Convective Star and WD=White 
Dwarf. Other paths also contribute but are less important than those shown. The Galactic 
disc age is chosen to be 9.5 Gyr 


Finall y, we would l ike to b ring the attention of the reader to the calcu - 
lations of IPuez et ^ IPuez and MathisI (I2OIOII : IPuez et ^ 

(and references therein) on the stability of non-force-free magnetic equilibria 
in radiative stellar regions (cores of Sun-like stars, envelopes of intermediate 
and high mass stars, and compact stars). Through the use of semi-analytic 
techniques they constructed and tested an axisymmetric non-force-free mag¬ 
netostatic equilibrium. Their numerical calculations recovered the instabilities 
which are characteristic to purely poloidal and toroidal magnetic fields and 
proved computationally, and for the first time, that only a mixed configura¬ 
tion is stable under all types of perturbations. Their important results can 
be applied to describe the magnetic equilibrium topology in stellar radiative 
regions (instead of choosing the initial field configuration arbitrarily) and to 
provide initial conditions for magneto-rotation al transport in state-of-the -art 
stellar evolution codes (e.g., see the models of iMaeder and Mevnetl 12014 for 
giants and supergiant stars). More general applications of these results to 
multi-dimensional magnetohydrodynamic computations are also envisaged. 


4 Neutron stars 

Magnetic fields of neutron stars at birth are estimated to be in the range 
~ 10^^ —10^^ G. However, measurements of the birth magnetic fields of neutron 
stars are always indirect. Sometimes the field or its dipole moment is measured 
today, e.g., from neutron star spin-down or Zeeman spectropolarimetry of 
suitable progenitors, and then extrapolated backwards or forwards in time 
respectively in the context of an evolutionary model. Otherwise, the birth field 
is inferred indirectly from some relic of the birth event, e.g., calorimetry of the 
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supernova remnant. If the internal magnetic field evolves slowly, so that its 
strength today approximately equals its strength at birth, then gravitational 
wave upper limits also provide constraints. We discuss each of these approaches 
briefly below. More information can be found in other chapters in this volume. 

Population synthesis models h ave a long history of being used to infer 


the birth fields of neutron stars (Gunn and Ostrikei 


1971)1 : Hartman et al 


19971 : Faucher-Gigucrc and Kas't^ 200til : Kiel et al.l 2008 1. Radio timing mea¬ 
surements of the spins, fields, and radio luminosities of the current neutron 
star population are combined with prescriptions for binary evolution, source 
kinematics in the Galactic gravitational potential, radio emission properties 
(e.g., beaming), supernova kicks, and observational selection effects to con¬ 
strain the spin and field distributions at birth. Population synthesis models 


lated and recycled 

, in large-sale radio multi-beam surveys (Manchester et al. 

2001 

: Morris et al. 

2002l:lKramer et al.ll2003l:lFaulkner et al.ll2004: Hobbs et al. 

2004 

), an updated Galactic electron density map (Gordes and Laziol2002fl. and 


proper motions from very-long-baseline interferometry. The results are that the 
birth spins are inferred to be normally distributed, with average birth period 
(Pq) = 0.30 s and standard deviation ap^ = 0.15 s, while the birth fields Bq 
are log- normally distributed, with (logHg) = 12.65 and fJ inp^ Bn = 0.55 {Bq in 
gauss) ( Faucher-Giguere and Kas^l2006l : Ixiel et al]l2008ll . The latter studies 


find no evidence for magnetic field decay over < 0.1 Gyr, assuming that the 
radio luminosity scales roughly as the square root of the spin-down luminos¬ 
ity. The no-decay conclusion sits in partial tension with the inference of field 
decay over ^ 10“^ yr in a different sub-po pulation of neutr on stars, e.g., the 


magnetars and th e thermal X-ray sources ( Pons et al.ll2012 ). Faucher-Giguere 


and Kaspi ( 200fill also find no evidence for bimodality in the distributions of 


Pq and space velocity. 

The energetics of the birth event, which are related to Bq and Pq, leave 
their imprint on the supernova remnant. From X-ray measurements of the rem¬ 
nant’s radius, temperature, and emis sion measure, one can infer the total blast 


energ y assuming Sedov expansion ( Reynolds and Ghevaliei 1981 : Reynold^ 


I 2 OO 8 II . The blast wave is powered by the core-collapse event, which depends 
weakly on Bq and Pq, and the relativistic wind emitted by the newly born 
pulsar, whose luminosity scales oc BqP^^. Drawing on X-ray Multi-Mirror 
Mission Newton (XMM-Newton) data, Vink and Kuiper ( 2006ll showed that 
the pulsar wind played an insignificant role powering the blast wave in two 
anomalous X-ray pulsars (AXPs) and one soft-gamma-ray repeater (SGR), im¬ 
plying Pq ^ 5 ms and hence no significant protoneutron star dynamo in these 
three objects at least. Steep density gradients in the interstellar medium can 
modify this con clusion, especially if the gradient correlates with the interstellar 
magnetic field ( Vigelius et akllioOTh . 

Gravitational wave upper limits and future detections furnish independent 
information on the internal magnetic field of a neutron star. At the time of 
writing, coherent searches for the Crab pulsar with the Laser Interferometer 
Gravitational Wave Observatory (LIGO) constrain its gravitational wave lu- 



























































16 


Ferrario, Melatos, Zrake 


minosity to be < 2% of the observ ed spin-down luminosity, thereby lim iting 
the internal field to Bq ^ 10^® G ( Abbott et al. 20101 : Aaa^^^ 2014 ). An 
analogous result has been obtained for the Vela pulsar (jAbadie et al.l 20 111 : 


Aasi et al.ll2014D . 


The gravitational wave strain scales as the square of the spin frequency, so 
the best constraints come from rapid rotators. Two classes of object are espe¬ 
cially interesting in this regard. First, the large deformation of a newly born, 
fast-spinning (Pq ^ 1 ms) magnetar caused by a super-strong internal mag¬ 
netic field (> 10^® G) radiates powerful gravitational waves, which should be 


detecta ble with Advanced-LlGO-class detectors up to the distance of the Virgo 
cluster dStella et ^ l2005t DaH’Osso and Stella l2007t IPall’Osso et al.l [2009t 


Mastrano et al.l2011 ; Melatos and Privmakl2014) . A future gravitational-wave 


detection of a millisecond magnetar would provide a direct measurement of Bq, 
if the distance is known (e.g., from a counterpart whose redshift is measured 
electromagnetically). Second, recycled millisecond pulsars have low external 
magnetic dipole moments, but their internal fields may be much stronger, 
i f the external dipo l e is r educed by accretion-driven diamagnetic screening 
( Pavne and MelatosI 12004 ). which leaves the birth field in the interior un¬ 
touched. Combining this scenario with rapid rotation and tight electromag¬ 
netic spin-down limits (ellipticity e < 10“® in some cases), recycled millisec¬ 
ond pulsars already yield some of the best constraint s on the internal Bp by 
any m ethod, ruling out Bq > 10^® G in some objects ( Mastrano and Melatod 
I2OI2I) . depending on the poloidal-toroidal flux ratio and whether the core is 
superconducting or not. 


The above conclusions change, if the internal field evolves significantly in 
ordina ry pulsars (as opposed to magnetars), cf. Faucher-Giguere and Kaspi 
( 20061) . Radio pulse evolution over long time-scales may be a signature of 
magnetic polar wandering or plate tec t onics, where t he field wanders inside 
the star ( Macv 1974 : Ruderman 199ll: iMelatod 2012 ). The recent discovery 
of changes in the flux ratio (0.1 per century) a nd component sep aration (0.6 
degrees per century) of the Crab’s radio pulses (ILvne et al.ll201.ll) is notewor¬ 
thy in this re gard, although other explanations like rad iative precessi on are 
also possible ( MelatosIl200(]t Barsukov et al.l[201^ I 2 OI 4 ]) . Stairs et al. ( 2000l) 
discovered switching between two discrete magnetospheric states (and hence 
two pulse shapes and spin-down rates) in PSR B1828-11. 


The central, enduring debate regarding the origin of neutron star mag¬ 
netic fields revolves, as it does in magnetic white dwarfs (see Sect.[3|), around 
whether the field is a fossilised relic of the progenitor’s field or is generated 
afresh by dynamo action in the protoneutron star in the first ^ 10 s after 
core collapse. Below, we summarise briefly the pr os an d cons of each scenario, 
drawing heavily on an excellent review by Spruit (|2009l) . The reader is directed 
to the latter reference for more discussion. 































































Magnetic Field Generation in Stars 


17 



log,„(M/Mj 


Fig. 5 The magnetic flux <pp for the most magnetic main sequence stars (squares), white 
dwarfs (circles), high fleld radio pulsars (triangles) and magnetars (stars) 


4.1 Origin of fields in neutron stars: the fossil hypothesis 


Although the surface magnetic fields in pre-main sequence stars, Ap/Bp, OB 
stars, magnetic white dwarfs and neutron stars cover a wide range of magnetic 
held strength values (from a few 10^ to 10^® Gauss), the magnetic huxes near 
the upper limit of the observed held ranges are all a few lO^^Gcm^. We 
show in Fig. [5] the magnetic huxes of the most magnetic compact stars and 
non-degenerate stars with radiative envelopes. This hnding has been used in 
support of the fossil held hypothesis for the origin of helds. 

The fossi l held hypothes is for the origin of h elds in neutron stars was hrst 
proposed bv IWoltiCT ( 1964h . Buderman (1972) also remarked that the simi¬ 
larity of magnetic huxes in magnetic white dwarfs and neutron stars could be 
explained through hux conservation during the evolution from main sequence 
to compact star stage. 

Ferrario and Wickramasinghd ( 2006b ) have investigated the effects of the 


fossil held hypothesis for the origin of magnetic helds in neutron stars by car¬ 
rying out population synthesis calculations for different assumptions on the 
distribution of the magnetic hux of massive (> SMq) main sequence stars 
and on the dependence of the initial birth period of neutron stars. They used 
the observed properties of the po pulation of isolated radio pulsars in the 1374 
MHz Parkes Multi-Beam Survey (IManchester et al.ll200ll) to constrain model 
parameters. These were then used to deduce the required magnetic properties 
of their progenitor stars. Their conclusion is that the fossil held hypothesis, 
which does not allow for magnetic hux loss in the post-main sequence evolu¬ 
tion, requires a very specihc distribution of magnetic helds for massive main 
sequence stars as shown in Fig.|6] for an assumed dipolar held structure. In 
this picture the held distribution in massive stars would be continuous and all 
massive stars would be magnetic. This distribution is predicted to have a peak 
near 46 G with low- and high-held wings covering a held range from 1 to 10“* G. 
About 8 per cent of main sequence stars would require to have helds in excess 
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Fig. 6 Predicted magnetic field distribution of maasive stars (8 — 45 according to the 

population synthesis modelling of iFerrario and WickramasinghS ll2006br) under the assump¬ 
tion of fossil fields and magnetic flux conservation 


of 10^ G and these would the progenitors of the highest-held neutron stars. 
This main sequence distribution of BO-0 type stars is an observable quantity 
and thus a prediction of the fossil held model. 


Further population synthes is calculations of the observed prope rties of 
magnetars were carried out by iFerrario and Wickramasinghel (l2008bl) under 
the assumption that magnetars originate from main sequence stars that are 
much more massive than those giving rise to normal radio pulsars. These 
studies were prompted by (i) the observation that most magnetars have been 
linked to progenit ors in the mass range ~ 20 — 45 Mq I Gaensler et al ] l2005l: 


Muno et al. 2006ll and (ii) the discovery of strong helds in massive 0 -type 


stars (see Sect. HD. Thus fossil magnetic huxes similar to those observed in 
the magnetars may already be present in stellar cores prior to their col¬ 
lapse to neutron star and would explain the presence of helds of up to ^ 
10^^ G in magnetars. According to Ferrario and Wickramasinghel ( 2008bll . the 
anomalous X-ray emission of magnetars would be caused by the decay of the 
toroidal held that does not contribut e to the spin down of the neutron star. 
Ferrario and Wickramasinghel ( 2008bll predict a number of active magnetars 


(see Fig. [ 21 ) that is consistent with the number of sources d etected by ROSAT 
and a Galactic birth rate compatible to that inferred by iKouveliotou et aT 


Since stars of spectral types BO to O are the progenitors of neutron stars, it 
remains to be seen whether there is a sufficiently large range of magnetic huxes 
in these stars to explain the entire range of magnetic huxes in neutron stars. 
So far, observations of magnetism in massive stars by the MiMeS, MAGORI 
and BOB collaboration s have led t o the d iscovery of a large number of new 
objects. The paper by Petit et id] ( 2013ll lists the physical, rotational and 
magnetic properties of about 60 highly magnetic massive OB stars. Most of 
the detected helds are in the > 10^ G regime with masses > 8 M 0 and radii 
> 5 R 0 . Under magnetic hux conservation, they would be the progenitors of 
the highest magnetic held neutron stars (> 10^^ G). 
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Fig. 7 Filled squares: observed magnetars; fille d circles: magnetars as derived from the mod¬ 
elling of iFerrario and WickramasinghS ll2008blh The solid line is an em pirically deter mined 
boundary (a “magnetar death line”) given by log(P) = 8.41og(P) — 20 llZhan i|2^ 


Problems with the fossil hypothesis 


Taken at face value, the fossil field arguments seem to be reasonable. Cer¬ 
tainly, there is no trouble accounting for the magnetar end of the distribution, 
since typical surface fields measured in Ap and Bp stars compress to give fields 
> 10^^ G in neutron stars. If anything, there is a problem at the lower end, 
where it is hard to produce neutron stars with B ~ 10^° —10^^ G without some 
other accretion-related physics. In fact, a main sequence star with a surface 
field B ~ 300 G and a radius R ^ 2Rq compresses to give B ~ 6 x 10^^ G. 
There are three problems with the fossil picture. First, the innermost ^ 1.4 Mq 
of the progenitor, which collapses to form a neutron star, occupies ^ 2% of 
the progenitor’s cross-subsectional area, reducing the maximum compressed 
field to B ~ 10^^ G. One can raise the limit somewhat, if the progenitor’s 
fie ld is centrally concentrated , but t he three-dimensional MHD simulations 
of Braithwaite and Nordlund ( 2006bll (which assume polytropic, spherically 
symmetric density and temperature profiles and neglect core convection) do 
not predict much central concentration. Second, t he magnetar birth rate is 
comparable t o the ordinary neutron star birth rate ( Keane and Krameill2008t 
Woodjl2nn8h . yet, as already mentioned in Sect. 14.11 there may be too few 


known progenitors with Bg P 10 kG to create magnetars in the numbers 
observed. Further more, the relatively low incidence of inter mediate-strength 
fields (10-100 G) ( Wade and the MiMeS Gollaboration 2014 ) may also repre¬ 
sent a problem for the fossil field hypothesis, unless future, more sensitive 
spectropolarimetric surveys of massive stars discover that intermediate fields 
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(e.g., see Fossati et al. 2014^ are indeed much more common than current ob¬ 
servations show. Third, the core of the progenitor couples magnetically to the 
giant envelope. The coupling is ‘frictional’ — the field lines are disrupted into 
a tangle by instabilities as they shear, behaving like a turbulent magnetic 
viscosity rather than ‘bungee cords’ — but it is strong nonetheless at high 
fields, scaling as Bg. Consequently, the core decelerates too much to explain 
the distribution of neutron star spin periods observed today. However, note 
that this is a problem for any neutron star with a magnetised progenitor, even 
if the neutron star’s field is generated ultimately in a dynamo; it is not an 
argument against fossil fields specifically. Indeed, a protoneutron star dynamo 
may never take root, if the progenitor core rotates too slowly. Core-envelope 
coupling has stimulated the suggestion that neut ron star spins are imparte d 
by the birth kick during supernova core collapse ( Spruit and Phinnevlll998h . 
Alignment is therefore expected between the spin axis and proper motion, for 
which some eviden ce exists in X-ray images of supernova remnants and pulsar 
radio polarimetry ( Wang et al.ll2006ll . 


4.2 Origin of fields in neutron stars: dynamo generated fields 


In this model, the seed field is amplified by a a dynamo in the proto-neutron 


star (Huderman and Sutherland 1973 ; Duncan and ThompsorJl992 ; Thompson and Diincan 


1993:lBonanno et al.ll2005). There are two main kinds of protoneutron star dy¬ 

namos: those driven by convection and those driven by differential rotation. 
Convection-driven dynamos were first analysed in the neutron star context 
by Thompson and Duncan (|l993h . A thermal luminosity L is transported 
outwards through the star by turbulent convection with characteristic eddy 
velocity V « (L/AirR^p)^^^, where p is the mean stellar density, and this gen¬ 
erates a magnetic field of strength Hpns ~ (47 rpF^)^/^, if the mechanical and 
magnetic stresses are in equipartition. For a typical photon luminosity, we ob¬ 
tain BpNS ~ 10® G (L/10®® erg s-i)i/®(A/10® cm)^/®. For a typical neutrino 
luminosity, Hpns is ~ 10® times higher (L ^ erg s“^). When the pro¬ 
toneutron star collapses from R = 10® cm to 10® cm, the field is compresses 
further by flux conservation to give B « IO^Hpns- 

A convective dynamo faces several obstacles to explain all neutron star 
magnetic fields, although it may certainly explain a subset. The most seri¬ 
ous obstacle is the assumption of equipartition. Numerical simulations reveal 
time and again t hat self-sustaining dynamos do not reach equipartition (e.g.. 
Cook et al.ll2003h because then the magnetic stresses would paradoxically self- 
quench the driving shear. In a variety of contexts, and for various numerical 
set-ups, the magnetic stresses satur ate at < 5% of the mechanical stresses 
( Cook et '^120031 : iBraithwaitel 2006ll . These theoretical results are consistent 
with observations of the Sun, whose dipole amounts to ~ 0.3% of equipartition 
(jCharbonneaul 120141) . A second obstacle arises when stratification quenches 
convection ^ 10 s after the dynamo begins. Does the field retreat or leave be¬ 
hind a permanent dipole moment? The danger here is that pockets (‘domains’) 
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of tilted polarity are left behind, if the convective eddies turn over faster than 
they are quenched, producin g a dipole moment lower than one inferred in mag- 
netars (e.g., from spin down iKonveliotou et al.lll998ri . even if the higher-order 
multipoles are substantial. The consequences for magnetic structure of forming 
convective pockets (which ‘insula te’) an d radiative pockets (which ‘conduct’) 
has been explored by Tout et al. ( 20041) . 

A second type of dynamo is driven by differential rotation. It is sometimes 
presumed that a sheared, magnetised fluid develops a toroidal field component 
which reacts back to shut off the shear completely. In practice, the back 
reaction can only ever be partial in a continuously driven system; if the fluid 
velocity gradient were to vanish momentarily, generation would cea se, and 
the g radient would eventually reassert itself in response to the driver ([Spruit 
I 2 OO 9 II . Microphysically, quenching is avoided because toroidal wind-up is lim- 
i ted to Brh 5, 10i3p by instabilities, where Bp is the poloidal field component 
(|SDrnitll2009l) . Hence it is necessary to grow B^ and Bp together i n a differen¬ 
tially driven dynamo. The d ynamo loop can be clos ed by Tayler ( Braithwaitel 
20061) or magneto-rotational ( Moiseenko et al. 20061) instabilities, for example. 
The growth window is open for ^ 10 s, before the star stratifies stably, and 
lasts for a time ~ R/va or {AS1)~^ in the latter scenarios respectively, where 
VA is the Alfven speed and Afi is the angular velocity shear. Differential ro¬ 
tation can also be driven by st ellar m ergers: the reader is referred to section 
and Wickramasinghe et al. ( 20141) for details. 

Braithwaitel ( 2006h showed that a shear-driven dynamo can even persist 


in a stably stratified star. Fluid is displaced parallel to equipotential surfaces 
in the Tayler instability. Hence, the condition for dynamo operation reduces 
to requiring that the maximum hydromagnetic mode wavelength allowed by 
buoyancy exceeds the minimum wavelength set by ohmic dissipation — a con- 
dition which is easily satisfied in a highly conducting neutron star. Braithwaite 
( 20061) showed explicitly through simulations that a Tayler dynamo is indeed 
self-sustaining in the presence of strong stratification and that it saturates 
with Bd, ^ 5 Bp a nd B'^/Stt Ri 10~^pV‘^, i.e., at 1% of equipartition. However, 
IZahn et ^ ( 2007 ) have strongly contested this result, finding that the Tayler 
instability cannot sustain a dynamo even at high magnetic Reynolds number. 

A mean-field variant of the shear-driven d ynam o was analysed in the neu¬ 
tron star context by Mastrano and Melatos ( 201l[) . In a mean-held dynamo, 
the turbulent electromotive force (v x B) in Faraday’s law is the sum of two 
terms: the ‘a’ term proportional to B, which arises from turbulent vortic- 
ity, and the so-called ‘Jl x J’ term, which is proportional to ? 7 i,(V x B)., 
and arises from off-diagonal entries in the resistivity tensor rjij (lMofFattlll978 : 
Biskamnlll997l) . A natural seat for a mean-held dynamo is a tachocline, w here 


the angular velocity jumps sharply ( Ott et al. 20061 : Burrows et al. 2007); in 
a protoneutron star, this may occur at the boundary between the convection 
and n eutron huger zones ( Bonanno et al.ll2003ll2005l) . Mastrano and Melatos 


( 2 OI 2 I) proved that a tachocline mean-held dynamo operates self-sustainingly 
in a protoneutron star, if is anisotropic and the jumps in ‘a’ (vorticity) and 
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angular velocity are coincident, and also (less realistically) if rjij is isotropic 
but the jumps are displaced. 


4-.2.1 Generation physics 

The ultimate fate of a dynamo-generated magnetic field depends on buoyancy 
— does the field float out of the star or remain anchored in the core? — and 
helicity — does the field maintain its large-scale integrity against the action 
of hydromagnetic instabilities, or is it shredded? 

Buoyancy exists in a neutron star because, in j3 eq uilibrium, the ratio of _ 

charg ed to neutral species densities increases with depth ( B.cisenegger and Goldreid] 
199211 . Magnetic flux tubes (‘light fluid’) rise to the surface if the Alfven speed 


satisfies va ^ HN, whe re H is the hydrostatic scale he i ght, and N is the 
Brunt-Vaisala frequency ( Reisenegger and Goldreich 1992 : Spruit 2009[l . The 
field rises most easily when the protoneutron star is young and neutrinos help 
to lower N to < 10 rad s~^. By contrast, in a normal neut ron star, the¬ 
ory predicts N ^ 10^rads“^ (jReisenegger and Goldreichlll992ll . In this way, 
magnetic loops inside the star penetrate the surface and add up to form a 
permanent dipole moment, as long as their footpoints remain anchored in the 
cooling crust or core and do not float away entirely. 

Assuming that the field does remain anchored, is it stable? Many analyses 
have shown that a purely poloidal magnetic field in a high l y conducting plasma 


nave snown tnat a purely poloidal magnetic neld in a nign l y conducting plasma 
is alw ays unstable (IMarkev and Tavlerllh?,? ; Wrightl[l973l : Flowers and RudermanI 
Il977ll . If we regard the star’s eastern and western hemispheres as two bar mag¬ 
nets, it is energetically favour able for the magnets to flip and counteralign 
destr oying the dipole moment ( Flowers and RudermanI[l97^ iMarchant et al. 


Honl); in practice, the instability proceeds by forming Rayleigh-Taylor ‘mush¬ 
rooms’ in the toro idal direction. Similarly, a pur ely toroidal magnetic field 
is always unstable ( Tavlei 1973 : Akgiin et al.ll20l5l . However, linked poloidal- 
toroidal field s appear to be stable under a range of conditions, b oth in analytic 
calculations (|Prendergastlll956t Tavleilll973 ; Akgiin et al.ll2013ll and numerical 
simulations dBraithwaite and Spruit 2004al : B^raithwaite and Nordhmdl 2006hl : 


Braithwaitel 2009ll . 


Linked poloidal-toroidal fields are stabilised, because their helicity is con¬ 
served approximately in a highly conducting plasma. Helicity is defined as 
H = fd^x A • B, where A is the magnetic vector potential; it is a gauge- 
dependent quant ity which descr ibes the number of times one flux tubes winds 
around another ( Biskamnl 1997 1. In a weakly resistive plasma, ff is conserved 
globally, even when iJ is scrambled locally on individual flux tubes. Hence, if 
a dynamo-generated field has H ^ 0 initially, the twist can never be undone 
completely, leading to H 7 ^ 0 in the final state. If magnetic energy is minimised 
during the fi eld’s evolution, it tends to a force-free configuration, cf. labor atory 
spheromaks ( Broderick and Naravan 2008t Mastrano and Melatosl[2008[) . 

Numerous papers have been wri tten recently reporting stability calcula- _ 

tions on poloidal-toroidal fields, e.g.. dBraithwaite and SpruitfeOOda ; Braithwa ite and Nordhmdl 


lions on puioiuai-ioruiuai neiyiSj_e;^j;jj|oramiwaiie_^u^_o;yruujizyynayDi;amiwaiie_^u^^Nurm 

2006b : Akgiin et al.ll2013 : Armaza et al.ll2013 : ICiolfi and Rezzollal2013l : lGourgouliatos et al 
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2 OI 3 I : Lander and ,Tonesl2012 : Mitchell et al.l20Ll Passamonti and LandeiJ2013 1 
A comprehensive review lies outside the scope of this chapter; we simply men¬ 
tion a few highlights. In a barotropic star, it has become apparent that sta¬ 
bility depends in subtle wa ys on the bounda r y conditions and the en tropy 
distribution within the star ( Ciolfi et al ] |201ll: [bander and Jonegir2012ll . In a 
non-barotropic star, magnetic field configurations with an equatorial ‘torus’ 
are ge nerica lly stable, when the poloidal-toroidal flux ratio is ‘modest’. Akgiin 
et al. ( 2 OI. 1 II obtained the stability condition 0.25 < B^/Bp < 0.5(7ad/7 — 
l)\Eg\/ (Bp/Sir), where 7ad and 7 are the adiabatic and non-barotropic specific 
heat ratios respectively, and \Eg\ is the total gravitational potential energy. 
Non-barotropic solutions also match self-consistently to an external dipole (or 
any arbitrary multipole), exte nding the power of fu t ure observational tests 
involving gravitational wav es ( Mastrano et al.l I2OIII: iMastrano and MelatosI 
2OI2I : [Mastrano et al.llioT^ . Finally, in a superconducting star, the stability 
problem changes character fundamentally for two reasons: the Lorentz force 
is much stronger (|Fb| oc HdB ^ where iLci is the critical field in 

a type II superconductor), and its vectorial structure is more complicated 
[Fb oc B X (V X Hci) -f pV(|B|i9|Hci|/i9p)]. The equilibrium magnetic struc- 
ture and its stability are then controlled by ibni / (B) at the crust-c ore boundary 
( Glamnedakis et ;il]l2012l : [bander et al.ll20ljlLandeil 201,11 2014ll . In fact, the 
internal field configuration now depends very strongly on field strength, sug¬ 
gesting signifi cant differenc es between the interior fields of pulsars and those 
of magnetars ( Landeil[2014ll . 


4-2.2 Superfluid turbulence 


Does magnetic field evolution in a neutron star conclude ^ 10 s after the 
protoneutron star is born and stably stratifies? The traditional consensus 
has been in the affirmative, except for t he slow ohmic a nd H a ll evolution 


descri b ed elsewhere in this volume and in iGennert et al.l (|2012fl : iPons et al 


( 2 OI 2 II : Viganb et al. ( 2012h . However, this picture may need modification in 
light of recent work suggesting that the neutron superfluid and charged compo¬ 
nents to which it is coupled are turbulent. In this scenario, internal magnetic 
activity may be ongoing, even though the magnetic flux thereby produced 
struggles to rise buoyantly to the surface. Hints of abrupt magnetospheric 
changes observed on ^ 1 yr time-scales (see Sect.|4]) — much longer than the 
dynamical time of the magnetosphere, but much shorter than the ohmic and 
Hall time-scales — are relevant in this context. 

In the standard picture, the vorticity field inside a neutron star is uniform 
and quasi-static. Macroscopically, this means that there is zero meridional 
circulation. Microscopically, it means that the superfluid is threaded by a rec¬ 
tilinear array of vortices, each carrying a quantum of circulation k = /i/(2m„), 
where h is Planck’s constant and m„ is the mass of the neutron. As the 
star brakes electromagnetically, an angular velocity lag AH builds up be¬ 
tween the crust and superfluid, whose value is set by the balance between 
the Magnus and pinning forces. One finds pnRAil « Fpi„(^coh'?pin)~^ and 
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hence AQ ^ 1 rad independent of Q and f?, where -Epin is the pinning 
energy, is the superfiuid coherence length, and ^pin is the pinning site sep¬ 
aration ( Warszawski and MelatosI 201,3li . Glitches reset AQ but only partially, 
because there is no observed correlat ion between glitch sizes and waiting times 


in most pulsars (IMelatos et al.ll2008ll . 


The lag Afi changes the flow structure completely. It is a well-known 
result of fluid mechanics that a shear flow in a sphere cannot be purely 
toroidal. The boundary conditions induce Ekman circulation, wherein inte¬ 
rior fluid is drawn into a viscous boundary layer, s pun down through contact 
with the crust, then recycled b ack into the interior ( Abnev and Epstein 1996t 


van Evsden and MelatosI 2010ll . This occurs with or without a rigid core. Fur¬ 


thermore, Ekman circulation is known to be unstable at high Reynolds num¬ 
bers (see Sect. E3) Re ^10^, forming unsteady, nonaxisymmetric flow patterns 
l ike herringbone wa ves (Re Ri 6 x 10^) and Taylor-Gortler vortices (Re ~ 10“^) 
I N akabavashill 19^ 1. In a neutron star, where the effective Reynolds number is 
much greater [Re > 10^ for e~-e~ shear viscosity modified by Landau damping 
by transverse plasmons; see Shternin and Yakovlev! ( 2008 )], wave modes res em- 
bling fully developed turbulence are expected ( Melatos and Peraltall2007ll . 

At the microscopic level, the vorticity field is disrupted by vortex-line insta¬ 
bilities. One example is the Donnelly-Glaberson instability, familiar fr om labo¬ 


ratory experiments with liquid helium and Bose-Einstein condensates ([Glaberson et al. 


19741) . If there is a superfluid counterflow directed along the vortices, Kelvin 


waves are excited and amplih ed by mutual friction to pro duce a reconnecting, 
self-sustaining vortex tangle ( Tsnbota et al.l 20031 l2013h . This process has a 


low threshold (counterflow Ri mms“^, easily e xceeded during Ekm an pump¬ 
ing) and a fast growth time (< one spin period) ( Peralta et al. 2006[). Recently , 
a related Kelvin-wave instability has been discovered by Link ( 2012a ; 2012bll , 
which does not rely on a counterflow and ari ses from imperfect p i nning , which 
creates a lag to amplify the Kelvin waves [cf. iGlamoedakis et al. ( 2009ll ]. This 
instability grows over ~ days, which is still fast. 

Turbulent, tangled vorticity is an important ingredient in the story of mag¬ 
netic field generation. The turbulent neutron condensate couples to charged 
components in the star through mutual f r iction and entrainme nt (IMendell 


llOOsI : IPrix et al 20021 : Haskell et al. 2009. 201^ Ghamell 2013ll. as well 


as 


through fluxoid-vortex interactions (jSrinivasan et al. 19MJahan-Miri 2010ll . 


Consequently, the charged components circulate too (jPeralta et al.ll2005 ). Sim¬ 
ulations that incorporate the magnetic field dynamics explicitly are required 
to understand the behaviour of the coupled fluids in detail. There is reason 
to be hopeful that such simulations will be conducted in the next few years, 
motivated by the prospect of gravitational wave experiments. 

Can superfluid turbulence be quenched by stratification, i.e., by the buoy¬ 
ancy force produced by the charge-neutral ratio increasing with depth (see 
Sect. l4.2TT|) ? Counterintuitively, perhaps, the situation is borderline. Strati- 
hed turbulence is characterised by two dimensionless quantities: the Froude 
number, Fr, which is related to the buoyancy force, and the Reynolds num¬ 
ber, Re ( Brethonwer et al. 2007ll . Even for Fr“^ > 10^, the turbulence is not 
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or 


quenched, if Re is sufficiently large. Turbulence persists for Re > Fr ^ 
equivalently R^(ASl)^ > where v is the kinematic viscosity (jlida et al 


2 OO 9 I : Laskv et al.ll201.‘i[h Neutron stars lie near the boundary defined by the 
above condition in the Fr“^-Re plane. In the regime Re > Fr“^ and Fr“^ 
> 10^, the flow is modified strongly away from Kolmogorov isotropy, even 
though it re mains circulatory; the Earth’s atmosphere and ocea ns also occupy 
this regime ( Brethonwer et akl 20071: Chung and Mathe^l2012h . 

Another possibility is that superfluid turbulence is quenched by magnetic 
stresses. Again, though, the case is unclear. At the micro scopic level , the im¬ 
perfect pinning instability which drives a vortex tan gle (Link! 120131) persists 
even when a magnetic field is present; in fact. Link’s ( 20131) analysis assumes 
magnetic locking of the viscous proton-electron plasma to the crust. Does the 
wound up by flux freezing erase the shear? More simulations are needed, 
but for now it seems unlikely: the magnetic back reaction satura tes at 5-. 1% 


of th e mechanical (shear ) stress in simulations attempted to date (ICook et al 


2003 : Braithwaitel boooh . and the Tayler instability self-organizes such that 
B • Vi? = 0, i.e., such that wind-up st ops, before quenches the pro¬ 
cess ( Dnez et al. 2006 : Rudiger et al.ll200^ . Magnetised Ekman pumping acts 
on the rapid time-scale ~ (i?i?/z;^)^/^i?“^ (where Va is th e Alfven speed ) 
in simple magnetic geometries, e.g., uniform magnetisation ( Easson 1979bl) . 
However, it is known to operate much more slowly in complicated g eometries 
and stratified conditions ( Goedbloed and Poedtsl[2004 : Melatosll2012 ). Finally, 
Easson showed that many — perhaps most — realistic magnetic topol ogies 
(e.g. , those with cl osed loops) cannot enforce corotation ( EassonI 'l979ali : see 
also Melatos ( 2012h . Physically, this occurs because the toroidal component of 
the magnetic field perturbation created by spin down satisfies a nonc onserva¬ 
tive eq uation of motion in a superconductor with Fermi liquid forces (lEasson 
1979ah . Again, large-scale numerical simulations are needed to determine if 


the class of magnetic topologies that permit corotation is restricted or generic. 


5 Status of MHD turbulence theory and computation 

The basic physics underlying stellar magnetic fields lies in the theory of mag¬ 
netised fluids, or magnetohydrodynamics (MHD). Of particular importance is 
the manifestation of magnetohydrodynamic turbulence because its presence 
underpins the operation of stellar dynamos, and controls diffusive magnetic 
flux transport in stars. This section is intended to give a high-level overview 
of recent progress in MHD turbulence theory and computation as it applies to 
stellar magnetism. 

5.1 Equations of motion 

The equations of compressible MHD can be written as 
dt 


V•(pu) 


(5) 
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i 9 u 

dB 


( 6 ) 

( 7 ) 


where in Equation [5] u is the fluid velocity, p is the mass density, and v is the 
kinematic viscosity. Equation[7]is Faraday’s law together with the phenomeno¬ 
logical Ohm’s law E = —u x B+77J. In general an equation for the transport of 
energy must also be supplied along with an equation of state relating the gas 
pressure to its temperature and density. However, if incompressibility is as¬ 
sumed, as it is in most analytical treatments of turbulent dynamo, the energy 
equation is superfluous. 


5.2 Dimensionless numbers 

If L and u represent characteristic length and velocity magnitudes of the sys¬ 
tem in question, then the ratio Re = wLjv is referred to as the Reynolds 
number. It represents the relative importance of inertial and viscous forces. It 
can also be seen as the ratio of the viscous time tvisc = jv to the advective 
time tadv = L/u. Flows having large Re are in general susceptible to turbu¬ 
lence as their dynamics are dominated by the non-linear coupling in Equation 
[HI Due to the enormous length scales involved in astrophysical flows, it is com¬ 
monly appropriate to drop the viscous term altogether and just work with the 
Euler equation. 

Analogously, one can define the magnetic Reynolds number Rm = uL /77 
which measures, correspondingly the relative importance of magnetic induc¬ 
tion and diffusion.lt is also the ratio tres/tadv where tres = is the damping 
time scale for magnetic fluctuations of scale L. The limiting case of infinite 
Rm corresponds to the case of a perfectly conducting fluid. This limit has the 
special property that the magnetic field is “frozen” into the fluid. Flux-freezing 
is a conservation law for the which the quantity B • dA. is a constant of the 
motion, when S represents a surface which is deformed along with the flow. 
The limit of vanishing viscosity and resistivity is referred to as ideal MHD. 

One can also define the magnetic Prandtl number Pm = Rm/Re which 
characterises the relative importance of viscous and resistive dissipation. In 
turbulent flows, it also corresponds to the ratio of length scales where the 
velocity and magnetic fields become smooth under the action of their respective 
dissipation mechanisms. Thus, for flows having large Pm, the velocity field is 
smoothed out by viscosity over scales where the magnetic field is still frozen 
into the fluid. 

As an example, let us estimate the magnetic Prandtl number for various 
parts of the sun. The resistively of a fully ionised plasma is given by IPriesd 
(I2OI4II 

me 

V = 


PoUeC'^Tf 


( 8 ) 
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where me is the electron mass, Ue is electron number density, and Te is the 
electron mean free path. In the solar photosphere this gives rj ^ 10"^ m^s“^. 
Thus vortices on t he solar surface of size ^ 500 km rotating with typical 
speed of lkms“^ ( Simon and Weiss! 1997^ have magnetic Reynolds number 
Rm ~ 5 X 10^. By comparison, the Reynolds number is of order 10^ making 
the magnetic Prandtl number of order 10^. A similar estimate for the solar con¬ 
vection zone gives Pm ~ 10“^. The magnetic Prandtl nu mber of diffuse astro- 

1015 


physic al plasmas is vastly larger, for examp le it is 


Maron and Cowlev 


( 200lll in the Galaxy and as large as 10^^ ( Schekochihin et al.l 2002b ) in the 
intergalactic medium. These values reflect their extraordinarily low collision- 
ality. 

Other dimensionless numbers of interest include the Lundquist number 
S = Lva/ii where va = is the speed of Alfven wave propagation. 

The sonic and Alfven Mach numbers are flow speed measured in units of their 
corresponding wave propagation speeds. The sonic Mach number Ai is also 
roughly the square-root of the ratio of kinetic to thermal energy densities, and 
similarly with the Alfven Mach number M.a except with the magnetic energy 
instead of the thermal energy. 


5.3 Conserved quantities 


Conserved quantities in MHD belong to one of three types ( BekensteiDlll987 ). 
The first kind, referred to as type A, is where the closed line integral of a vec¬ 
tor field is a constant of the motion, where the path of circulation is carried 
along with the local flow. The canonical example is that of Kelvin’s circulation 
theorem, which states that in a perfect fluid the circulation F = u • ds is 
a constant of the motion, where the closed path C is advected and deformed 
along with the flow. Kelvi n’s cir c ulatio n theorem has been generalised to rela¬ 
tivistic h ydrodynamics by Taud (|l95^. and to non-relativistic and relativistic 
MHD by iBekenstein and OronI ( 200Clf) . Another example is the Alfven law of 
magnetic flux conservation, which states that A • ds is a constant of the 

motion. This law can also be stated as the constancy of magnetic flux through 
the surface whose boundary is the closed contour C. It formally applies in ideal 
MHD, where the viscosity and resistivity are vanishing. This law, commonly 
referred to as “flux freezing”, is generally considered a good approximation 
in astrophysical plasmas, where the microphysical diffusion time of magnetic 
fields is far longer than other dynamical times. However, for turbulent flows 
in which the velocity field stays chaotic down to scales below those of interest 
(referred to as “rough” velocity fields) turbulent diffusivity effectively violates 
the perfect fluid assumptions. The violation of Kelvin’s circulation theorem in 
high R eynolds number hydrodynamic flows has been analyse d by IChen et al 
( 200fill , and similarly for the Alfven theorem by lEvink ( 201 ill . 


The second type of conservation law (type B) is where a local fluid quan¬ 
tity is a constant of the motion along streamlines. Bernoulli’s theorem is the 
canonical example of type B conservation laws. 
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Type C conservation laws apply to the global conservation of scalar or 
vector quantities whose transport is strictly local. These are like Gauss’s law 
in that they relate the conserved quantity’s time rate of change in a finite 
volume to its flux through the volume’s bounding surface. They can thus be 
stated in either differential or integral form, the latter of which forms the 
basis for modern numerical methods for solving the hydrodynamic and MHD 
equations. The prototypical type C conservation law is the continuity equation 
p + V • pu = 0. In the absence of viscous and resistive effects, momentum and 
magnetic flux are also conserved in this way, where pu is replaced with a 
generic flux function for that conserved quantity. Similarly, in the absence of 
heat transport by conduction, the evolution of total energy can also be written 
as a continuity equation. 


5.4 Consequences of helicity conservation 


A particularly important conserved quantity is the magnetic helicity Hm = 
/ (PxA ■ B. Hm is a global topological characterisation of the magnetic field, 
which indicates the degree to which magnetic field loops are mutually linked. 
It is analogous to the kinetic helicity Hk = f (Pxu-ui (where a; = V x u is the 
vorticity) which represents the degree to which the vortex lines are mutually 
linked. In a perfect fluid, both helicity measures are a constant of the motion. 

Magnetic helicity is important because it tends to establish magnetic fields 
well above the turbulence integral scal e. This consequence of magnetic helicity 
conservation was first investigated by Frisch et al. ( I975h . where it was pre¬ 
dicted that the magnetic energy would peak around the scale Hm/Eb (where 
Eb is the total magnetic energy), even if the injection of magnetic energy took 
place at much smaller scales. This process is referred to as inverse cascading 
because the magnetic energy cascades from small scales toward large scales, 
unlike the direct cascade of energy observed in hydrodynamic turbulence. It 
is the underlying mechanism for the formation of large-scale dynamos which 
may be responsible for establishing stellar magnetic fields. 


5.5 Inverse cascade of magnetic energy 


Magnetohydrodynamic turbulence exhibits an unusual phenomenon where, at 
least in its freely decaying state, magnetic field correlations arise over dis¬ 
tances much longer than the original eddy size. This inverse cascading effect 
is reminiscent of what is seen in two-dimensional (but not three-dimensional) 
hydrodynamic turbulence. Alt hough this effect m ay be related to the con¬ 


servation of magnetic helicity ( Frisch et al. I975h . until recently it was not 


appreciated that it happens even when the field is initially helicity-free. That 
is, when it is topologically equivalent to the uniform (or zero) magnetic field. 

Inverse cascading of ma gnetic energy was first seen in numerical simulations 
bv lMenegiizzi et aLl dQSlTl . It has since been confirmed by many others using 
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direc t numerical simulations IcThristensson et al. 20W; Baneriee and Jedamzikl 


2004 : Kahniashvili et al. 2010 : iTevzadze et al. 2012), quasi-analytical t ech¬ 
nique s such as the “Eddy damped quasi normal mode” approxima tion (|Son 
I 999 I) and also turbulence “shell models” ( Kalelkar and Panditll2004l) . Such in¬ 
verse cascading is a process of magnetic field self-assembly, whereby magnetic 
energy provided by random motions at very small scales, may over time estab¬ 
lish dynamically important magnetic fields over scales that are much larger. 
For example, the sun’s convective cells, which set the turbulence integral scale, 
are far smaller than the solar diameter. The ordered magnetic fields that ex¬ 
tend over the whole surface of the sun (and many other stars) may come from 
such a self-assembly process. Dynamos capable of establishing magnetic fields 
on scales far larger than the turbulence integral scale are generically referred 
to as large scale dynamos. 

The literature to date is still conflicted on whether non-zero magnetic he- 
licity is a nec essar y condition for in verse cascading to occur. It was shown by 
OlesenI (1997) and Shiromizul ( 19981) that inverse cascading could be expected 
even for non-helical magnetic field configurations, as a consequence of rescaling 
symmetries native to the Navier-Stokes equations. But no inverse cascading 
was seen in numerical studies based on EDQNM th eory (TS onl 199 91) or direct 
numerical simulations with re latively low resolution ( Christensson et al ] [20nt 
Baneriee and Jedamzikl 20041) . 

Very recently, three numerical studies have appeared which confirm the 
existe nce of inverse cascading i n non-helical three-dimensional MHD turbu¬ 
lence. Brandenburg et al. ( 20141) carried out high resolution direct numerical 
simulations of compressible MHD turbulence with magnetic energy injected 
at very small scales. The results co nfirm the gro wth of magnetic energy spec¬ 
trum PM{k, t) at low wavenumbers. [Zrakel ( 2014 ) found similar results in high 
resolution simulations of freely decaying relativistic MHD turbulence. In the 
latter study, the evolution of PM{k,t) was parametrised around a self-similar 
ansatz which may be used to evaluate the time required for magnetic fields 
to be established over s cales arbitrarily larger than the energy injection scale. 
Berera and LinkmannI ( 2014h also found inverse cascades in decaying MHD 
turbulence. That study utilised a novel technique, whereby the system was 
evolved hydrodynamically with a passive magnetic field. The resulting energy 
spectrum also exhibits growth of magnetic field fluctuations at large scales, 
suggesting that inverse cascading of passively advected vector fields is a generic 
property of hydrodynamic turbulence. 


5.6 Small-scale turbulent dynamo 

In general, the term “dynamo” refers to the conversion of kinetic energy into 
magnetic energy. In magnetohydrodynamic theory, a number of distinct pro¬ 
cesses can facilitate this conversion. The only thing they all have in com¬ 
mon is the presence of turbulence. Many excellent reviews have been writ¬ 
ten on this the connection between MHD turbulence and dynamo processes: 
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Brandenburg and Siibramanian (2005); Lazarian and Cho (2005); Knlsrnd and Zweibell 
Tobias et alJ ( ^llh . Here we will present a very brief overview of a 


process known as small-scale turbulent dynamo. 

Small-scale turbulent dynamo operates in highly conducting turbulent flu¬ 
ids, even when the turbulence is isotropic and non-helical. Conceptually it is 
quite simple. Magnetic field lines carried along by the turbulent flow become 
increasingly chaotic and distorted over time, leading to enhancement of their 
energy. When the field is very weak, i.e. va u, its back-reaction on the 
fluid can be ignored, and the field evolves as a passively advected vector field 
in otherwise hydrodynamic turbulent flow. This approxim ation is referred t o 


as kinematic small-scale dynamo, and was first studied 


( 1968 1 . 

(l2nnili . 


but more modern analytic treatments exist such as that of IVincenz 
The original formulation, known as the Kraichnan-Kazantsev dynamo starts 
with the assumption of an arbitrary spectrum of velocity fluctuations that 
are delta-correlated in time. Despite its simplicity, this model makes accurate 
predictions regarding the time rate of change and spectral energy distribution 
PM{k,t) of the magnetic field, where k is the inverse length scale. In partic¬ 
ular, solutions for PM{k,t) that are exponentially growing in time exist when 
the microphysical diffusion coefficient rj is small relative to the turbulent dif- 
fusivity, a condition which is met when the Reynolds number of the flow is 
sufficiently large. 

The Kraichnan-Kazantsev dynamo predicts for the unstable growing solu¬ 
tions, (1) that PM{k,t) oc and peaks at the resistive cutoff, and (2) that 
the magnetic energy exponentiates at the fastest eddy turnover time of the tur¬ 
bulent hydrodynamic cascade. The first numerical experime nts to confirm the 
basic features of the Kraichnan-Kazantsev dynamo were bv iMeneguzzi et al 


( 198lh . At the time, direct numerical simulations were exceedingly expensive, 
and those simulations were limited to grid sizes of 32^ and 64^. Since then, 
many groups have continued this investigation utilising high performance par¬ 
allel computing architectures, and modern numerical schemes for obtaining 
solutions to the MHD equation. 

Although the original work by Kazantsev predicted that the magnetic 
Reynolds number Rm had to be “sufficiently large” for small scale dynamo 
action to occur, the exact value of the critical R m needed to be found t hrough 
numerical experiments. For the case of Pm = 1, iHaugen et al.l (|2004bll found 
that the critical magnetic Reynolds number Rmcrit for dynamo action was 
about 35. For smaller Pm, it was found that Rmcrit increases as Pm decreases. 
This may be partially understood as having to do with the flow properties in 
the resistive interval. When the Pm is very large, the resistive scale lies deep 
below the viscous cutoff, so that magnetic energy peaks where the flow is very 
smooth. However, when Pm is smaller than one, the resistive scale moves into 
the turbulent inertial range of the fl ow, and the velocity field is rough at the 
resistive scale. iHaugen et al.l (|2004cl) foun d that Rmcrit, inc r eased to 70 when 
the flow is supersonic, even when Pm = 1. iFederrath et al.l ( 201lh found that 
for compressible flows the dynamo growth rate depends upon how the turbu¬ 
lence is forced; vertically driven flows produced more efficient dynamo than 
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did di latational forcing. It was argued analytically bv iBoldvrev and Cattane^ 
(|2004ll that in the small Pm regime, the critical magnetic Reynolds number 
should be a strong function of the roughness exponent a of the velocity field, 
oc such that smaller a (a rougher field) requires a larger Rm to produce 
small-scale kinematic dynamo. 

A central question around small-scale dynamo is the nature of its sat¬ 
uration. The exponentially growing solution of the Kazentsev model holds 
only as long as the dynamo is in the kinematic regime, so that magnetic 
back-reaction on the flow is negligible. Once the field attains a strength high 
enough to compete with the fluid inertia, the kinematic assumption breaks 
down, and some sort of saturation behaviour is expected. This saturation 
is a highly nonline a r prob lem, and thus very difficult to study analytically. 
Schekochihin et al. ( 2002all argued that in the non-linear stage, at least for 
very large Pm, magnetic energy would remain concentrated in the sub-viscous 
range, moving toward the inertial range on a resistive rather than dynamical 
time. Were this to be the case, equipartition of the turbulent kinetic and mag¬ 
netic energy densities at the largest scales would require an asymptotically 
long time to be established. However, this scenario has not been observed 
in the numerical literature. Instead, the emerging consensus (other differences 

asidel dCho and Vishniac 2000: Brandenburg et al.ll200^: Haugen et al.l2004a: 


Schleicher et al.l20Rll:ISchekochihin et al.l2004l : Maron and Cowlevl200lt Zrake and MacFadven 


2011 , 2013h is that magnetic energy shifts into the inertial range regardless of 
the magnetic Prandtl number. Once there, the field attains coherence at in¬ 
creasing scales until finally reaching the energy-containing scale of turbulence, 
thus establishing scale-by-scale equipartition. Both the kinetic and magnetic 
energy follow a Kolmogorov spectrum, with the magnetic energy per unit 
wavelength exceeding that of the kinetic energy by a factor of around 2. 

The time-scale for the flow to fully esta blish scale-by-scale equipartition is 
still not completely clear. It was argued bv iBeresnvakI (2012) that non-linear 
small-scale dynamo exhibits universality in the sense that a constant fraction 
Ce oi the work done by turbulent pumping accumulates in magnetic energy. In 
this scenario, scale-by-scale equipartition should be atta ined after a time \ICe- 
Carefu l modelling of the approach to saturation by Zrake and MacFadvenI 
(201 ,lll finds that indeed scale-by-scale equipartition is reached after a number 
of dynamical times that is independent of the Reynolds number. 


6 Conclusions 

In this chapter we have reviewed progress made on the origin of magnetic fields 
in intermediate and massive main sequence stars and in compact stars. 

Observations of magnetic fields in white dwarfs may indicate that the dis¬ 
tribution is bi-modal, exhibiting a paucity of objects between 10^ — 10® G, 
although further sensitive spectropolarimetric surveys conducted on 8 m class 
telescopes are needed in order to confirm this finding. Should this bimodality 
be confirmed, then two different channels for field formation and evolution 
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could be at play. In particular, one could identify weakly magnetic main se¬ 
quence stars such as Vega and Sirius as the progenitors of weakly magnetic 
white dwarfs belonging to the low field component. This would lend some 
support to the fossil field hypothesis for the origin of fields in compact stars. 

However, the viability of the fossil field hypothesis has been recently ques¬ 
tioned because of the absence of magnetic white dwarfs paired with non¬ 
degenerate stars. These observations mirror the dearth of magnetic non-degenerate 
A, B and O stars in short-period binary systems. It is not clear whether there 
is any connection or analogy to be made between fields on the main sequence 
and compact star phases, but a merger scenario for held formation has been 
suggested to overcome some of the problems raised by recent observations. 

Here we stress that the pairing properties of magnetic white dwarfs do 
not necessarily invalidate the viability of the fossil scenario since more than 
one formation channel could be at work. However, the complete lack, rather 
than just a dearth, of magnetic white dwarfs paired with non-degenerate com¬ 
panions in non-interacting systems, raises strong doubts on whether the fossil 
route can be the main channel for the formation of magnetic white dwarfs. If 
magnetic white dwarfs were the descendants of magnetic main sequence stars, 
then magnetic white dwarfs should be commonly found in binary systems to 
reflect the incidence of binarity observed in (magnetic) stars of F to late B 
spectral types. However, there is not one single example of a magnetic white 
dwarf paired with a non-degenerate companion in a non-interacting binary. 

In order to explain this curious finding, it has been proposed that fields 
in isolated highly magnetic white dwarfs could be generated during common 
envelope evolution when the two components of a binary system merge. 

The observational situation surrounding birth magnetic fields of neutron 
stars remains frustratingly uncertain, because what measurements exist are 
indirect. Population synthesis models point, on balance, to moderate birth 
fields ~ G. Zeeman spectropolarimetry has revealed strongly mag¬ 

netised (i.e., a few thousand gauss) progenitors across a range of spectral 
classes, but even so it seems that the numbers may be insufficient to account 
for the magnetar birth rate under the fossil field scenario. At least for magne- 
tars, therefore, a proto-neutron star dynamo is favoured, driven by differential 
rotation an d the Tayler o r magneto-rotational instabilities rather than thermal 
convection ( Spruit 2009ll . Emergence (through buoyancy), stability (through 
helicity), and ongoing activity (through superfluid turbulence) play roles in de¬ 
termining the strength of the (observable) surface field. Future gravitational 
wave observations offer the best prospects for direct internal field measure¬ 
ments. Until these observations become available we can only say that, as for 
magnetic white dwarfs, it is possible that both fossil and dynamo processes 
are at work in different neutron stars. 

We have also reviewed the theory of magnetohydrodynamic turbulence 
because its presence underpins the operation of stellar dynamos and controls 
diffusive magnetic flux transport in stars. But MHD turbulence poses a great 
theoretical challenge due to its inherent complexities, and many features of 
its operation are yet to be satisfactorily explained. Advances in numerical 
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algorithms and computing resources have yielded considerable gains in the 
empirical grounding of various phenomenological pictures. 
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